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Abstract 



o 

^ 1 The interior of a neutron star is likely to be predominantly a mixture of superfluid neutrons and 

. superconducting protons. This results in the quantisation of the star's magnetic field into an array 

■ of thin fluxtubes, producing a macroscopic force very different from the Lorentz force of normal 

\£) [ matter. We show that in an axisymmetric superconducting equilibrium the behaviour of a magnetic 

field is governed by a single differential equation. Solving this, we present the first self-consistent 
superconducting neutron star equilibria with poloidal and mixed poloidal-toroidal fields, also giving 
the first quantitative results for the corresponding magnetically-induced distortions to the star. The 
Mh poloidal component is dominant in all our configurations. We suggest that the transition from normal 

to superconducting matter in a young neutron star may cause a large-scale field rearrangement. 

Introduction 



O 

1 There is now compelling evidence that the bulk of a neutron star is composed of superfluid neutrons and 

superconducting protons. It has long been thought that neutron stars would be cold enough to contain 
these states of matter [TJ [2J [3], based on expectations from the theory of terrestrial superconductivity [JJ. 
A long-standing piece of evidence in favour of superfluidity is the phenomenon of pulsar glitches - 
■ sudden events in which the star spins up slightly [5]. The best explanation for the larger glitches is as a 

£NJ , transfer of angular momentum from the more rapidly rotating superfluid to the crust. In addition, recent 

observations of the rapidly cooling young neutron star in the Cassiopeia A supernova remnant are well 
explained by the onset of neutron superfluidity together with proton superconductivity in the core [5] [7] , 
suggesting that these components will be present in a neutron star below a critical temperature of around 
10 9 K. The most highly magnetised neutron stars, the magnctars, have high observed temperatures which 
might suggest their transition to superfluidity is delayed. This is due to crustal heating however, whilst 
the thermally decoupled magnetar core is expected to cool rapidly below the critical temperature jS]. 
Thus all observed neutron stars are likely to contain superfluid and superconducting components. 
In terrestrial superconductors the Meissncr effect expels any magnetic field below some critical strength, 
but in neutron stars the expulsion timescale is extremely long [2], so the magnetic field exists in a 
metastable state. Nonetheless it is affected: in a large part of the star the protons are expected to form 
a type-II superconductor, meaning that the field will be quantised into fluxtubes surrounded by umnag- 
netised matter. On the macroscopic scale this changes the nature of the magnetic force from the Lorentz 
force of normal matter to a fluxtube tension force [9] [TUl E] ■ The nature of the innermost part of the 
core may exhibit type-I superconductivitv|12 [ [T3j . with large regions of alternating normal and supercon- 
ducting matter; the effect of this on the global magnetic field is unknown at present. Alternatively, the 
star's hadronic matter could give way to an inner core of 'colour superconducting' quark matter [T3I [To*] . 
A neutron star's magnetic field can play a variety of important roles. It affects the temperature and 
rotational evolution of the star — and rotation is the key observational feature used to determine the 
star's age. Understanding the interior fields of apparently different classes of neutron star could help 
unify them into a single canonical stellar model [T7J [16] . Magnetic field effects could explain the differing 
nature of glitches in pulsars and magnetars and post-glitch recover v[T5). In magnetars, the magnetic 
field provides the energy that powers their giant flares and is important for understanding their observed 
quasi-periodic oscillations [HI [20]. Finally, a magnetic field induces a distortion which will generally 
not be aligned with the star's rotation axis. This system will therefore produce gravitational waves [21] . 
perhaps at an amplitude great enough for future detection. 

Motivated by the above reasons, there has been a great deal of recent work on neutron star magnetic 
fields, but almost all of it based on models assuming normally-conducting matter. This may be partly 
because this case is more familiar, thanks to the large body of literature on non-degenerate stars [22]. For 



neutron stars superconductivity is a key missing ingredient, whose inclusion is an essential step towards 
more realistic models. 

In this letter we try to lay some foundations for the modelling of superconducting neutron stars. We 
study equilibrium configurations, motivated by the observation that neutron star magnetic fields appear 
to be long-lived, evolving only on long timescales. We show that in axisymmetry the magnetic field of 
a superconducting star is governed by a single differential equation. This is analogous to the 'Grad- 
Shafranov' equation for normal matter[531 [23], but more complicated: it involves terms related to the 
local field strength. We describe a method of solution for the equation and present the results. Other 
than the special case of a purely toroidal field[551 [55], these are the first self-consistent solutions for 
magnetic fields in a superconducting star (see, however, the simplified poloidal-field model constructed 
in |27j). We compare with the corresponding normal- matter results and discuss the implications for the 
internal magnetic fields of neutron stars. 



Axisymmetric magnetic stars 

We model a neutron star as a three-fluid system, with superfluid neutrons, electrons and type-II su- 
perconducting protons. The electrons have negligible inertia however, and their chemical potential may 
be incorporated into that of the protons. Therefore we can reduce to a two-fluid system of equations, 
denoting neutron quantities with a subscript 'n' and the combined proton and electron quantities with 
a 'p'. We choose a rather idealised equation of state, a two-fluid analogue of a polytrope [28], so that 
the chemical potential /i x of each species (x € {n, p}) is a function of the corresponding mass density /3 X : 
Mx = Ax(Px)- Working in cylindrical polar coordinates (w,<f>,z), the two Euler equations for our model 
may be written: 

where $ denotes gravitational potential, £7 X rotation rate and F x magnetic force. Although we only 
consider non-rotating models here, the following derivations and numerics follow through for cases with 
corotating neutrons and protons, f2 n = f2 p . In general one would expect an entrainment effect, leading 
to coupling between neutrons and protons and to an effective magnetic force on the neutrons. We neglect 
this for simplicity, so that F n = 0, and write the proton force as ^mag. The particle species are therefore 
only coupled through gravity: 

V 2 $ = 4 7 rG( /9n +p p ). (2) 
Now, taking the curl of the proton Euler equation, we see that there exists a scalar M such that 

^ = VM, (3) 
Pp 

regardless of whether the protons are normal or superconducting. Another universal result is that B 
must be divergence-free; using this together with the assumption of axisymmetry allows us to write 

B = — Vu x e + B e , (4) 

ZD 

which defines the streamfunction u. Note that B ■ Vu = 0; field lines are parallel to constant-u contours. 



Normal matter 

In this familiar case $ mag is the Lorentz force, given by 

S r mag = ^(VxB)xB. (5) 

It can be shown that the toroidal field component is governed by a function fx (u) = wB^ (S3] [21] ■ Using 
this and the general results of the last section, one may derive the Grad-Shafranov equation: 

d 2 u I du 8 2 u 2 dM df N 

A * U = :5 _ 2 A h ~a~2 = -^TTppW JN-7-- (6) 

ovj z w dm oz z du du 



This governs the form of a poloidal or mixed poloidal-toroidal field in axisymmetric equilibrium. The 
single-fluid version of this has been the basis of numerous studies of magnetic stars. 



Superconducting matter 



By averaging the contribution of the fluxtubcs in a typc-II superconductor, one arrives at a macroscopic 
expression for the magnetic force [9l [TUl [TTj : 
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dH 

B ~df~ 



(7) 



where H cl = i/ c iB is the 'microscopic' critical field, B the smooth-averaged 'macroscopic' field and 
B = B/B the unit tangent vector to the magnetic field. In the absence of entrainment H c \ = h c p p to 
a good approximation, where h c is a constant Using this relation and equations Q and (|7|), we get 
an expression for the magnetic force with various poloidal terms (see |26j ) and one toroidal term, which 
must be zero by axisymmetry: 

(&nag)* = Vm X \J{p p wB^) = 0. (8) 

This gives us a dichotomy: satisfying ((SJ leads cither to a mixed poloidal-toroidal field or a purely toroidal 
field. For the latter case, we satisfy © by taking Vu = 0; this has been discussed in earlier work [2"5ll2l)] . 
For the former, we require that Vu and \7(p p vuB^) be parallel, which leads to 

PpwB^ = f(u) (9) 
for some function /. In the special case f{u) = the field is purely poloidal. 

One key step in the derivation of the Grad-Shafranov equation is showing that M = M(u). In the 
superconducting case this is no longer true. We can, however, define a related function y which is a 
function of u: 

, \ 4ttM , , 

y(u) = ——+B. (10) 
h c 

Using the functions y(u) and f(u), together with the poloidal magnetic- force terms from |26j . we arrive 
at a single differential equation governing the magnetic field: 

vn.y« 2 dy 2 df , . 

A*u = — w p p u- — n /— , 11) 

11 au du 

where II = B/p p . The above equation is the equivalent of the Grad-Shafranov equation when the 
protons form a type-II superconductor instead of being a normal fluid. The result for a single-fluid 
superconducting star may be obtained by replacing p p with p, the total density, and is valid for barotropic 
equations of state. The most significant differences from the normal case are the presence of the II factors 
and the fact that the magnetic force no longer appears explicitly through the function M. 



Superconducting core and normal 'crust' 

The superfluid and superconducting matter of a neutron star's core does not extend to the stellar surface; 
instead the star has an elastic crust of normal matter. It is also numerically difficult to solve equation 
(fTTj) matched directly to a potential field (V x B = 0) exterior. For these reasons we choose a canonical 
neutron star model with a core of superfluid neutrons and type-II superconducting protons, matched to 
a single-fluid relaxed 'crust' composed solely of normal protons; this in turn is matched to a vacuum 
exterior with potential field at the stellar surface 

The complicated physics at a neutron star's crust-core boundary may well include the presence of a 
current sheet and a discontinuity in the magnetic field. Quantifying these effects is beyond the scope of 
this paper, however, so for this first study we assume force balance and continuity of the magnetic field 
at the boundary. In particular, we define the crust-core boundary as an isopyenic contour p p = p p c (at 
a near-constant radius of r = 0.9i?»), requiring that the magnetic- force scalar function M be continuous 
and that H c i — > B there. This may be done by defining the superconducting functions y(u),f(u) in 
terms of their normal-matter counterparts Mn(u), /at(u): 

y{u) = h cP c p c + 4irM N (u)/h c , (12) 



/(«) = f N (u)/h c . (13) 
This then allows for smooth matching of the right-hand sides of the two governing equations © and (1111) . 
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Figure 1: Structure of a poloidal magnetic field in a superconducting star. On the left we show the 
magnitude of the field and on the right its direction (i.e. the field lines). The stellar surface i?„ is 
indicated with the solid arc at a dimensionless radius of unity, whilst the dashed line at 0.9-R* shows the 
location of the crust-core boundary. 

Numerics 

The governing equation for superconducting equilibria is harder to solve than the Grad-Shafranov 
equation ([6]) . The latter is itself unusual in having the argument u appear on both left- and right-hand 
sides, but the superconducting version also has the quantity II, which implicitly involves derivatives of u: 



With direct solution seeming infeasible, this problem is suited to a numerical iterative method, where both 
left- and right-hand sides are gradually updated until the scheme converges and produces a consistent 
solution for u. We use an adapted self-consistent field method [IS [26], employing an underrelaxation 
step for the solution of (fTTj) . 

Using the virial theorem [30] we confirm that our results are indeed equilibrium solutions, with a relative 
error of the order 10~ 5 . This error decreases accordingly with increasing resolution. As another check, 
we have constructed equilibria with increasing field strength (and hence distortion), confirming that the 
induced distortion scales in the expected manner: linearly in H cl B [3_T][!|]. Details of these tests and the 
numerical method will be presented in a later paper. 



We consider a stratified stellar model with a composition gradient, i.e. the proton fraction p p / p varies 
within the star. Various forms for the magnetic functions are permissible. We choose /at(m) = i(w~Wint) < ' 
when u > u- mt and /jv( u ) = otherwise, where u; nt is the largest contour of u (i.e. field line) that closes 
within the star; this avoids having an exterior current. We also take Mn(u) = ku 2 unless otherwise stated. 
Here a, ( and k are constants related to the strength of the magnetic field components. The figures are 
presented in dimensionless units, since the structure of the magnetic field is essentially independent of its 
strength. 

We begin by looking at a purely poloidal field, in figure [TJ This is broadly similar to the corresponding 
field in a normally-conducting star, with the field strength attaining a maximum J5 ma x deep in the star 
and vanishing in the centre of the closed field-line region. The main difference is that the superconducting 
star shown has a far larger 'weak-field' region than a model with normal protons; more specifically, a 
larger volume of the superconducting star has a field strength B < 0.1i? max . This effect is even more 
pronounced for models with Mn(u) = ku. 

In figure [5] we show a typical mixed poloidal-toroidal field configuration. For stars with normal pro- 
tons the toroidal component fills the weak-field region shown in figure [I] producing a 'twisted-torus' 
configuration [26]. This is partially true here too — but at the centre of the closed field line region, where 
the poloidal field vanishes, the toroidal component vanishes too. The resulting toroidal-field geometry 
is hence tubular (in three dimensions). The mixed- field configuration in figure [21 like all those we have 
found, is dominated by the poloidal component; the toroidal component only contributes 0.7% of the 
total magnetic energy. This may be related to the fact that equation (|II[) has a purcly-poloidal limit but 
no purely-toroidal one. 
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Results 




Figure 2: A mixed poloidal-toroidal magnetic field in a superconducting star. We show the poloidal 
field lines in black and the magnitude of the toroidal-field component with the colour scale. The stellar 
surface and crust-core boundary are indicated with the solid and dashed circles. Note that the toroidal 
component has a tubular structure, vanishing in the centre of the closed field line region. 



Table 1: The ratio B/B s of average internal field to polar field and the prefactor k e of the elliptic- 
ity formula (|T5|). for different poloidal-field configurations (specifically, different relations between the 
streamfunction u and the magnetic- force function Afjv(u)). 
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1.5 
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2.5 


k e [lQ- 8 } 


2.5 


3.4 


3.9 


3.5 



Next we look at the magnetically-induced ellipticity e = (<5oq — Qpoic)/Qeq, where Q cq , Q po \c represent the 
components of the star's quadrupole moment along the equator and pole. We rescale to a typical 1.4-solar 
mass neutron star with a radius of 10 km and assume a purely poloidal field. The mass of the normal-fluid 
crust is small, so the ellipticity scaling is well approximated by that of a purely superconducting star: 

« = 3.4X10-^^, (15) 



,10 12 G J V10 16 G, 

where B s denotes the surface field strength at the pole. We adopt a central critical field of H c i(0) = 10 16 
G, using the approximate formula given in |llj . For comparison, in the same stratified two- fluid model 
but with normal protons we find (using the code described in [26]): 

A simple way to approximate the ellipticity of a superconducting star is to take a result for normal matter 
and scale it up by a factor H c i/B, taking H c i = 10 15 G. Comparing our ellipticity formulae (fl~5j) and (fT6|) . 
both from self-consistent calculations, we see that for a given B s this approach would underestimate the 
star's distortion by around 30%. Choosing Mn(u) as a higher power of u increases the ratio of average 
internal field B to B s and hence the ellipticity at a fixed value of B s . This is summarised in table [T] 



Discussion 



We have described a method to solve for the magnetic field in a neutron star with typc-II superconducting 
protons. The magnetic force is more complicated in this case than for normal matter, but for an axisym- 
mctric equilibrium we show that it may be simplified to a single differential equation in the streamfunction 
and the local field strength, in analogy with the Grad-Shafranov equation of normal matter. We solve 
this using an iterative scheme, presenting the first self-consistent models of a superconducting neutron 
star (other than the special case of a purely toroidal field). 

Perhaps the most notable difference from normal-matter models is the generic appearance of a region in 
the star where the field strength vanishes. This may have repercussions for the dynamics of a neutron 



star; in normal matter it is known that such a region leads to an instability [32] [53]. Whether such an 
instability occurs in superconducting matter is an interesting open question. 

The interior field of a neutron star cannot necessarily be inferred from the observed exterior dipole field. 
If many poloidal field lines close within the star, or if there is a strong toroidal component, the internal 
field could be much stronger than expected from outside — a 'hidden' energy reservoir for the star. In 
our models the average interior field strength is 1.5 — 2.7 times that at the polar surface. The contribution 
of the toroidal component appears to be generically small, however; in our equilibria it accounts for less 
than ~ 1% of the magnetic energy. This is different from recent simulations for main-sequence stars, 
which found stable equilibria with large toroidal components |34j . 

A leading theory for the origin of neutron star magnetism (for magnetars in particular) is that dynamo 
action in the young star generates a strong, dominantly toroidal field [35] — in contrast to our equilibria 
for mature neutron stars. This may cast doubt on the validity of our equilibria or the dynamo scenario. 
Alternatively, both could be reasonable — then, as a hot young neutron star cools to a multifluid state 
with superconducting protons, its strong toroidal field would no longer be in equilibrium. It would have 
to undergo large-scale rearrangement to a poloidal-dominated equilibrium, a potentially violent transition 
which could be observable. 

The effect of superconductivity on neutron star magnetic fields has been, to date, neglected by the vast 
majority of studies. This letter demonstrates that in the simplest equilibrium models it may be accounted 
for using similar techniques as for normal-matter stars. Two key issues which future work should address 
are the presence of a magnetic force on the neutrons and the physics at the crust-core boundary. Beyond 
this, superconductivity will surely also play an important role in the evolution and dynamics of neutron 
stars. 
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